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Debris discs and comet populations around Sun-like stars: 
the Solar System in context 
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ABSTRACT 



Numerous nearby FGK dwarfs possess discs of debris generated by collisions 
among comets. Here we fit the levels of dusty excess observed by Spitzer at 70 |J.m 
and show that they form a rather smooth distribution. Taking into account the tran- 
sition of the dust removal process from coUisional to Poynting- Robertson drag, all the 
stars may be empirically fitted by a single population with many low-excess mem- 
bers. Within this ensemble, the Kuiper Belt is inferred to be such a low-dust example, 
among the last 10 % of stars, with a small cometary population. Analogue systems 
hosting gas giant planets and a modest comet belt should occur for only a few per cent 
of Sun-like stars, and so terrestrial planets with a comparable cometary impact rate to 
the Earth may be uncommon. The nearest such analogue system presently known is 
HD 154345 at 18 pc, but accounting for survey completeness, a closer example should 
lie at around 10 pc. 

Key words: planetary systems - circumstellar matter - infrared: stars 



1 INTRODUCTION 

Impacts pose a hazard to life on Earth, especially in the 
case of an event enerReti c enouRh to destro y much of the 
surface crust and oceans (jZahnle et al.ll200it . The most de- 
structive impacts are associated with comets rather than as- 
teroids, as the former h it at hi^h speeds w hen infalling from 
the Kuiper Belt region (|jeffers et al.ll200ll ). but the present- 
day cometary impact rate is low, largely because most of 
the p rimordial bodies have been dispersed (|Morbidelli et al] 
This is thought to have occurred in the first Gyr of 
the Sun's life, when Jupiter and Saturn crossed a mean mo- 
tion resonance, destabilising the gas giants so that Saturn 
had close encounters with Uranus and Neptune. The expan- 
sion of their orbits perturbed the primordial Kuiper Belt 
jGomes et al.ll2005l ). producing the Late Heavy Bombard- 
ment of the Earth at around 700 Myr, after which the comet 
population was greatly depleted and the impact frequency 
has been much lower. However, the presence of the gas giants 
is still significant today, as dynamic al interactions can bring 
comets into the inner Solar System - iHorner fc Jone"^ (|2008l ) 
have recently shown that the impact rate on the Earth would 
vary significantly if the mass of Jupiter were different. 

Thus, in extrasolar systems we may expect that the 
number of comets and the architecture of the giant-planet 
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system will strongly affect the impact rate on any terrestrial 
planets present. Direct evidence of impacts in the inner sys- 
tems comes from detections of dust grains at temperatures 
of a few hundred Kelvin, representing break-up debris from 
parent planetesimals. Such detections via a mid-infrared ex- 
cess signal above the stellar photosphere are rare, but this 
does not imply that inner system planetesimals and plane- 
tary imp acts are spars e, as the debris lifetime at a few AU 
is short llWvattllioOsh . Less directly, we can examine the 
far-infrared debris signatures to assess the numbers of bod- 
ies colliding within the cool outer comet belts, and consider 
the influence of gas giants in perturbing comets inwards to 
where they threaten terrestrial planets. 

Here we use the results of recent volume-limit e d sur - 
vey s of Sun- like star s with Spitzer by iTrilling et al] (|2008l 'l 
and lBeichman et al.l (|2006l ) to study the dustiness of Solar- 
analogue systems, via the predominantly 70 \im detections 
of excess above the stell ar phot o spher es. Various analytical 
models, as reviewed bv IWvattI (|2008l ). predict far-infrared 
dust luminosities as functions of initial size and mass of the 
planetesimal disc and the age of the host star. Thus the ex- 
cesses detected by Spitzer can be related to the ensemble 
of populations of comets per star. We can then place the 
Sun's Kuiper Belt, for which there are measurements of the 
numbers and distribution of comets and dust particles, in 
the context of similar stars. For example, other systems are 
known that are hundreds of times dustier than the Solar Sys- 
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Figure 1. Sun-like stars in order of decreasing 70 \aa excess (y-axis). The rank on the x-axis is the star's number (e.g. 1 for the most 
dusty) divided by the number of stars surveyed. Upper limits are not plotted. The diamonds show contiguous debris detections while the 
four points plotted with circles have intervening systems with only upper limits in E-jq ; the spacing of these sparse detections is derived 
as described in section 3.1. The black curve is a power-law fit for stars of Ejq down to 0.5, excluding the four sparse detections of lower 
excess. The grey line is an extrapolation to lower excesses including dust removal theory (section 3.4), not a fit to the four circle symbols. 
The dashed line shows the 2.5(t upper limit to the net E-jq for stars where detections beyond x = 0.3 could be made (section 3.2). The 
circle symbol (far right) is the position of the Sun in the ensemble, derived from its maximum estimated dust level (section 4.1) and 
plotted at its corresponding lowest rank on the grey curve. The bars indicate that the Sun may be less dusty and so of higher rank. 



tern; dustier at greater age than the Sun; or have much larger 
debris belts. Viewed externally as a debris disc, the Kuiper 
Belt would be faint, a result which motivates this study to 
answer the question of whether the planet and planetesimal 
population around the Sun is unusual, and whether this has 
any role in impacts on the Earth and the evolution of life. 



2 DATA 

Two major unbiased surveys of nearby Solar analogues 
have been made with Spitzer, with debris detection rates 
of «15 % at 70 (im, for fractional excesses of a few 
tenths or more above the photosphere (iTrilUng et al.ll2008l : 
iBeichman et all 120061 '). The FGK targets in these two sam- 
ples are complementary and extend out to ~ 30 pc from 
the Sun; a further survey (|Kim et al.l |2007| ) will com- 
plete much of thi s volume. We neglect here results from 
the FEPS project jmUenbrand et al.ll2008l : ICarpenter et all 
I2OO9I : I Meyer et al.l 20081 ) where stars were selected by age. 



and so there is a trend of lower sensitivity for younger and 
thus rarer and more distant stars. The published nearby- 
volume samples also have some biases that are inherited 
from the goals of the original proposals, with more luminous 
stars, hosts of giant planets and single stars over-represented 
compared to their true proportions among Su n-like stars. 

Only 70 M-ni data of high q uality from iTrilUng et al.l 

l|2008l ) and IBeichman et all l|200d ) are discussed here, trac- 
ing cool dust at tens to hundreds of AU from the host 
star. The photospheric signal was not detectable with 3- 



sigma confidence in all the stars, so there is a potential bias 
where faint stars rise above the detection threshold of to- 
tal flux only in cases of large excess. To avoid this, we se- 
lected only the objects with F,/SF > 3, where -F* repre- 
sents the photosphere and 5F is the noise in the flux mea- 
surement. This selection includes 176 stars from the two 
combined surveys, after also eliminating a number of M- 
dwarfs from Beichman ct al. (2006). The excesses are then 
defined as E70 = Fdust/F, = [Ftotai - F^)/ F,. All further 
analysi s uses the fluxes and errors listed by I Trilling et al.l 
(200^): [Beichman et al] (|2006l ) along with their other tabu- 
lated quantities, such as estimates for stellar ages. 

The survey papers identify debris systems as those with 
-E70 at 3(7 confidence, i.e. X70 = Fdust/SF > 3, here com- 
prising 27 stars. Here we add 6 more candidate systems of 
slightly lower significance. A tabulation of all 176 E-jo val- 
ues, in order of X70 , shows that negative outliers extend down 
to — 2.5cr, and this distribution is quite symmetric, with 70 
stars of ^ xto ^ +2.5 versus 73 of —2.5 ^ X70 < 0. Fur- 
ther, this X70 distribution is close to a Gaussian centred on a 
null value, with mean and sigma of -0.03 and 1.2 respectively. 
We thus adopt here as candidate debris systems all those 
with > -1-2. 5(T significance, and estimate that the chance of 
a false positive is < 1.5 % per star, given that < 1 more ex- 
treme outliers appear among 73 negative results. With this 
new selection criterion, there are 33 systems identified here 
as debris host^ , or 19 % of the stars observed. 



1 Those of 2.5-2.9(7 confidence are HD 39091, 55575, 69830, 
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Figure 2. Debris properties of the sample, in logarithmic age 
bins. The grey bars show the minimum-to-maximum E^o within 
each bin, and the points (and Poisson error bars) indicate the 
percentage of stars with debris. The number of systems detected 
(observed) in each age bin from left to right is 3 (12), 9 (45) and 
19 (99), with 18 stars having no tabulated ages. 



3 ANALYSIS 

3.1 Detected population 

Figure 1 plots the 70 \im excesses of the stars where debris 
is identified, ordered by a fractional 'rank'. The rank is for- 
mulated as the number of the star in order of decreasing 
E70, divided by the total number of stars (176). It is equiv- 
alent to the fraction of stars with equal or higher E70 than 
the current plotted point. However, beyond a rank of 0.1, 
not all the stars were observed to sufficient depth to detect 
the corresponding E70 of approximately 0.8 and below, if 
such excesses are present. This tendency is in fact implicit in 
the initial selection, as SEro = SF/Fi, and so where F,/SF 
approaches 3, only E7o > 0.83 can be detected with 2.5a 
confidence. In this regime, there are four debris detections 
(circles in Figure 1) that are interleaved with upper limits 
of similar Ero, but the ranks of the non-detections can not 
be determined as the true debris levels are unknown. How- 
ever, plotting these four points as contiguous ranks causes 
an artificial downturn in the trend of excess, because simi- 
lar systems should occur in this region but were not actually 
observed deeply enough to be detected. Thus for these four 
points, estimated ranks are plotted, by increasing the step on 
the X-axis from 1/176, to 1/176 multiplied by a scale factor 
of (number of stars still to be counted) divided by (number 
of stars where the next £70 could be detected). These stars 
with En ~ 0.4 were numbers 30 to 33 in order of decreasing 
dustiness, and their original ranks of 0.17 to 0.19 (30/176 to 
33/176) have been increased to 0.21 to 0.30. 

The excesses observed range from 17 down to 0.4, al- 
though only about one in five stars was observed deeply 
enough to detect debris at this lower bound. The minimum 
detectable excess for Spitzer, where generally the dust and 
star lie within the same telescope beam, is set by uncer- 
tainty in the photospheric predictions plus errors in the 
shorter-wavelength data used for extrapola tion. Dispersions 
in 24 ^.m data of around 5 % are observed (|Beichman et al.l 



114613, 173667 and 222237; HD 69830 is a confirmed deb ris sys- 
tem with shorter-wavelength excesses llTriUing et al.ll2008l'l . 
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Figure 3. Detected debris systems, showing stellar properties. 
Small, medium and large symbols correspond to the early, mid 
and late age-bins respectively of Figure 2, and light and dark 
grey shadings pick out the most and least luminous stars (early- 
F and K types respectively). Points within circles are systems 
hosting giant planets (from the Extrasolar Planets Encyclopedia) 
and those within squares have multiple stars (from the Catalog 
of Components of Double and Multiple Stars). 

l2006h and these are likely to be the minimum inherited er- 
rors when attempting to detect low levels of debris at 70 fi.m. 

There is no strong trend of level or incidence of debris 
with particular stellar properties. For e xample, there is littl e 
trend with age, a s alrea d y noted by e.g.iTrilling et aP l|2008l l ; 
iBeichman etal] (|2006l ): ICreaves fc Wvatj (|2003h . Figure 2 
shows the sample divided into three broad bins each span- 
ning a factor of fa 4 in age (to minimise errors from poorly- 
dated objects). The detection rate is essentially flat with 
age over the entire main sequence, while maximum dusti- 
ness declines mildly with age - a similar decline has bee n 
seen for 70 ^.m excesses of A- type stars (ISu et al.l |2006| V 
These tren ds agree with recent evolutionary models as sum- 
marised bv lWvat"^ (|2008l ) , where each system's dustiness de- 
clines slowly with time as the parent population of colliding 
bodies is ground away. If the initial planetesimal discs have 
a wide range of radii and masses, then there is a wide range 
of excess values at any one age, and the detection rate is 
rather constant with time if the upper env elope of maximum 
dust fiux is well above the detection limit l|Greaves &: Wva"t3 
l2003h . 

Further, there are only weak trends associated with 
other stellar properties, such as l uminosity, binari t y and 
the presence of giant planets |Beichman et al] l2006l : 
iBrvden et a"l]|2006l : rivming et al.ll2007l ). Figure 3 replots the 
debris hosts of Figure 1 with symbols denoting these prop- 
erties, and no obvious trends are visible. In further analysis, 
all the stars are therefore treated as a single population. 

3.2 Trend in debris level 

The ranked plot (Figure 1) shows that the excesses form a 
rather smoothly declining distribution, dominated by lower 
values of Ero- The steep slope can explain how the rate 
of far-infrared debris detection among Sun-like stars has 
risen with improvements in sensitivity with successive in- 
struments, namely the IRAS, ISO and Spitzer missions. For 
example, a typical G-dwarf photosphe re at 15 pc has a flu x 
of approximately 30 mjy at 60 ii.m (|Habing et akl 1200 ll ) , 
while the IRAS all-sky survey could detect signals of around 
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for the mean value and standard error on the mear0. This 
is not a significant detection of n et excess, so the susses - 
tion of a typical excess of 0.05 by|A umann fc Good! (|1990D 
(at 60 \xm) can not be confirmed. However, it does allow all 
systems to be slightly dusty, or some to have modest excess 
while others are dust-free, as well as the possibility that all 
are dust-less. Among these 24 targets there may in fact be a 
slight skew towards more positive E-jq measurements (Fig- 
ure 4), but a Kolmogorov-Smirnov test finds a probability of 
only 54 % that this positive tail are from a different (i.e. de- 
bris disc) population. Therefore there is as yet no conclusive 
result, but this does not rule out many low-dust systems. 
As a small excess is characteristic of the Solar System, at- 
tempting to infer the underlying population is useful to put 
the Sun in context. 



Figure 4. Cumulative distributions of measured E-jq (X axis) for 
24 upper-limit systems discussed in the text. These have small 
i?70 errors of ^ 0.15 and so the distribution is a guide to whether 
a net positive excess may be present. Positive values (solid line) 
are compared to negative values that have been multiplied by -1 
(dashed line). 



200 mjy (|Wvattl l2005') and the corresponding hmit for ISO 
was around 90 mJy ( Dccin et al...2000 ). Thus IRAS and ISO 
could detect iJgo ~ 6 and 3 respectively. Of the stars plot- 
ted in Figure 1, we would then expect about 3 bright exam- 
ples to be discovered by IRAS, and HD 139664, 48682 and 
109085 were actually found. Subsequently, ISO could have 
discovered ~ 6 less-dusty systems, while the sparse surveys 
actually undertaken added two discoveriefl, HD 30495 and 
17925 

The large sensitivity improvement with Spitzer has 
more than tripled the discoveries of debris systems among 
these nearby stars, suggesting that more fainter exam- 
ples remain to be discovered. However, debris levels in the 
presentl y undetected popu lation are inferred to be small on 
average. lAumann' fc Good (|1990) suggested that a net ex- 
cess of ^ 0.05 might be present, from IRAS 60 \im data for 
G-stars within 25 pc. The Spitzer 70 \ira data now explore 
this regime lying at the right of Figure 1, although individ- 
ual excesses are still not detected. For example, there are 24 
stars with Spitzer observations deep enough that excesses 
ranging from 3 % to 37 % could have been detected, if such 
dust levels were actually present in these systems. (Here we 
neglect the uncertainties due to calibration and just consider 
signals of 2.5 times the noise level.) These stars are part of 
the population at star-rank > 0.3 in Figure 1, and none of 
them has a significant measurement of excess, but the data 
are deep enough that a tendency to mildly positive values 
can be searched for. In particular, stacking the data can show 
if there is evidence for an average non-zero dustiness that is 
not apparent in each individual noisy observation. Averag- 
ing the measured £70 for these 24 stars we find 0.03 ± 0.02 



^ Five other fainter debris systems in Figure 1 were also in 
fact candidates from IRAS/ISO, where fainter excesses could be 
probed for various reasons, such as stellar proximity, or hotter 
photospheres or warmer dust than average. 



3.3 Single population hypothesis 

We first test the simplest possible hypothesis, that all Solar- 
type stars have debris discs drawn from a single distribu- 
tion. A suitable function is sought that fits the excesses, 
plus information in the low net excess for systems that are 
not detected. The single-population hypothesis is unlikely 
to explain all systems completely, but provides a testable 
empirical approach. 

Figure 1 shows (black line) a best fit to the debris de- 
tections with E70 = 0.5 to 17, given by 

£70 = 0.09a;~'"** (1) 

where x is the star rank, and the correlation coefficient is 
0.98. The single-population hypothesis is thus a rather good 
empirical fit to the systems detected, in spite of the wide 
range of stellar ages, spectral types, etc. The regime at E-jq < 
0.5 was excluded as a different dust behaviour is predicted 
(see below), and includes only 4 detected systems. 

If the Equation 1 fit continued for all ranks up to 1, 
then the last star would have E70 of 0.09, and the mean 
in the undetected population at a; > 0.3 would be 0.16. If 
this extrapolation applied, more detections should have been 
made, given a 3-sigma threshold of ~ 0.15 set by uncertain- 
ties ~ 0.05 in the photospheric levels as discussed above. 
Also, the true mean is unlikely to be so high, as among 
the 24 stars with the smallest excess-errors, there are only 
four objects with E70 measurementts even as high as ~ 0.25 
(2.5ct upper limits ~ 0.4). This suggests a drop in dustiness 
among the majority population of stars. 

3.4 Dust generation 

The hypothesis above assumes one distribution for the 
amount of dust per star, but given some range of initial cir- 
cumstellar disc properties, it is more likely that the quan- 
tity of comets per star forms the underlying distribution. 
In this case, the amount of dust generated needs to take 

^ The net value from all stars without debris detections is simi- 
lar at —0.01 ± 0.02 but is more influenced by unphysical negative 
excesses, down to E70 = —0.7. These may indicate uncertainties 
in the photospheric models, with the two lowest values corre- 
sponding to stars at the extremes of the spectral type range, for 
example. 
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int o account the dus t remo val processes, and as discussed 
by iDominik fc Jiooj '). this depends on the number 

of coUiding bodies. In massive discs, the dominant destruc- 
tive mechanism is colUsional and Ndust oc Ncomets, while 
in diffuse discs where the Poynting-Robertson (light drag) 
effect removes dust, Ndust oc Ncomets- This distinction be- 
tween high- and low-mass discs, in the coUisional and PR- 
drag dust removal regimes respectively, is also discussed by 
iMever et all (l2007l ). 



Wvatd dSoi) has shown that the discs readily de- 



tectable in the far-infrared are in the high-mass regime, but 
to extrapolate here to the whole population of stars, the 
dust generated in the low-mass cases must also be calcu- 
lated. This transition occurs where the fractional luminosity 
of the debris is 



fcrit 



t/L, = [2 X 10\r/dr){r/M,y^^]- 



(2) 



for disc radius and width r, dr in AU and stellar mass M* 
in solar masses. It is obtained from the limit where PR- 
drag and collision lifetimes are equal, for gra ins where th e 
radiative force is half the gravitational force (Wvatt 2005). 
Using 11 resolved images of debris discs aro und Sun-like 
stars observed in the optical and submillimetre l|Wvattl2008l : 
iGreaves et aP l2009l ). we calculated their values of fcrit- 
These lie in the range 1 to 10 xlO~® with an average of 
5 X 10~®. Th e corresponding 70 [i m ratio Ecrit can then 
be calculated (jSeichman et al.ll2006t ). and for dust tempera- 
tures of 40-80 K and stellar effective temperature similar to 
the Sun, Ecru is on average 0.5 at 70 \im. Hence, if all Sun- 
like stars possess debris discs of similar scales to the resolved 
system^ then there should be a downturn in dustiness be- 
low E70 ~ 0.5, even if there is a single smooth distribution 
of number of comets per star. 

The dotted line in Figure 1 shows this extrapolation be- 
yond £70 = 0.5, assuming that the remaining systems are 
less dusty according to i? oc Ndust oc Ncomets but continu- 
ing the same power-law distribution for Ncomets(x), i.e. of a 
form like Equation 1. In this case, the predicted mean E-jo 
is reduced from 0.16 to 0.06, in better agreement with the 
estimate from the data of 0.03 ± 0.02 for all the stars at 
X > 0.3. Thus, the single-population hypothesis taking into 
account the different regimes of dust removal can match all 
the present observational results. 



are similar at 60 and 70 jxm but the photosphere is fainter 
by a factor of 0.5±0.1 in the Spitzer data, from 6 G-dwarfs in 
common and with no debris. Hence, we adopt E-jq = 1.5 for 
r Ceti, while for the parameters of this system E^rit ~ 0.85, 
and so -ErCcti ~ Ecrit X 1.8. Since r Ceti ha s an estimated 
1.2 M ffi of colliding bodies of radii r ^ 25 km (jCreaves et al.l 
|2004) . under the assumption of a coUisional cascade with 
incremental population dN(r) oc r~^'^, then Ecrit would 
correspond to 0.7 M®. However, models (|Lohne et al.i 2008. 
e.g.) consider colliders up to r = 75 km, so including this 
extra contribution raises the critical mass to 1.3 M^. Thus, 
with the assumption that the dust in all systems comes from 
a coUisional cascade fed by 75km objects, the general scaling 



S70/O.5 = {Mr^7Skm/l.3 M®)" 



(3) 



where 0.5 is the generic value of Ecrit, and q is 1 for E > 
Ecrit and 2 for sub-critical discs. 

The mass in bodies of r ^ 75 km per star is then pre- 
dicted to range from 60 down to 0.2 Earth masses, with the 
detected debris discs hosting > 1 M® of comets. The shape 
of the population is the same as the excesses fitted in Fig- 
ure 1, i.e. Mcomets approximately inversely proportional to 
X, so the most massive belts are rare. The probability that a 
star's comets exceed a given mass is found by inverting the 
rank-mass expression to give 



P{> Mc, 



0.26/ A'Icomcts 1 



(4) 



applicable to the mass in bodies with r 5C 75 km. The high- 
est mass appears realistic, for example it is comparable to 
the 2 0-30 Mffi of colliders i nferred around the A-star Fomal- 
haut J Wvatt fc pentll2002ll. which hosts a bright debris disc 
with Evn > 20 (|Stapelfeldt et al.l 120041 '). The exact masses 
would be modified if t he comets do not f ollow a coUisional 
cascade - for example. iHahn fc Mahot"ral (|2005l ) fit the size- 
distribution for the Kuiper Belt as a double power-law with 
many small bodies but a depletion at sizes greater than 
65 km. Adopting this distribution for the exo-comet belts 
would increase the fraction of the total mass that is in the 
collider population by a factor of about 2.5. 



3.5 Comet populations 

The number of comets per star can next be estimated from 
the dust excess, albeit via several scalings so the the final 
values are only approximate. To connect comets and ex- 
cesses, we consider the r Ceti debris disc, which is one of 
the nearby re s olved examples and has E^o = 0.75 ± 0.15 
iHabing et al.l I2OOII ). These ISO observations were made 
with a slightly wider and shorter wavelength-centred fil- 
ter than the 70 \un. Spitzer/MIPS data, so we estimate a 
correction using stars detected by bot h ISO and Spitzer 
ijHabing et al ] |200ll : iTriUing et al.l [20oi ). The dust signals 

* Only 3 of the 11 resolved systems are in our present sample, 
and there could be a bias to larger systems where the aim is to 
make resolved images, although the resolved examples were not 
chosen for observation on such a priori expectations. 



3.6 Models 

Although the fit used here is simply empirical, analytical 
models for evolving debris disc populations in fact pro- 
duce similar trends. The models o f IWvatt et all (120071 ) have 
been applied to FGK stars surrounded by a range of plan- 
etesimal discs. These discs have masses drawn from a log- 
normal function and sizes from a power-law distribution, as 
in the previous mo del for debris discs of A-type host stars 
(| Wvatt et al.|[2007l ). The model systems are sampled as if 
observed at random ages and then ranked as in Figure 1. 
The general behaviour of a steep decline from a few very 
dusty systems down to a long tail of low-debris stars is eas- 
ily reproduced, adopting planetesimal populations of a few 
Earth masses and discs extending out to a few tens of AU. 
Kains et al. (in prep.) will present the final model parame- 
ters fitting both 70 and 24 \im. debris statistics among the 
Spitzer-surveyed FGK stars. 
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4 THE SOLAR SYSTEM 

4.1 Kuiper Belt properties 

We now compare the 70 iJ.m excess of the Kuiper Belt (if 
it were viewed from outside) to the extrasolar values. The 
Solar System's far-infrare d dust flux is unfortu nately not 
well established, although iLandgraf et"ai] (|2002l ) conflrmed 
that Kuiper Belt dust does exist, via impacts onto the Pi- 
oneer spacecraft as they entered the outer Solar System. 
However the emission from this dust from an Earth view- 
point is_a^_small_2erturbation to the much brighter Zodiacal 
flux.|A umann fc G ood ( 1990,) used IRAS scans to estimate 
that the peak Kuiper Belt flux is 3 x 10~* W m~^ sr~^ at 
60 \ira (in a 40 iJ.m bandpas s), and the main belt on the 
sky is effectively ^ 10° widfEl jMorbidelli et al.ll2003l '). Sum- 
ming over this area and dividing by the bandpass to give 
a per-frequency-interval unit then gives a total flux of ap- 
proximately 1 MJy (±0.5 MJy from the scatter in different 
scans). Subsequent COBE data with a similar 60 ^im pass- 
band favour an uppe r limit at the lower end of this range 
(jSackman et al.lll995l ). We assume that the 70 (xm dust flux 
will be similar, while the Solar photospheric emission would 
be 100 MJy if seen at the approximately 45 AU distance of 
the Kuiper Bel10. 

The Sun's excess at 70 |J.m is therefore « 0.01 ± 0.005 
from IRAS, or possibly less from the COBE limit. These 
values of excess for the Kuiper Belt are remarkably small. 
Allowing a 2(t uncertainty so that £70 can be as high as 0.02 
would imply that x ^ 0.9. Hence, the Sun is one of the least 
dusty systems (Figure 1), or possibly not a member of the 
single population being tested in our hypothesis. 

The total mass in Kuiper Belt bodies has been es- 
timated at a few hundredths of an Earth-mass, from 
deep surveys finding objects as small as 25 km diame- 
ter (B ernstein et al. 2004). The two-power-law model of 
iHahn Mahotral (j2005j) can be used to estimate the to- 
tal belt mass, adopting their generic albedo of 0.04. The 
Kuiper Belt would then sum to ~ 0.08 in the r ^ 75 km 
regime (or somewhat smaller in a coUisional cascade, or with 
more recently-adopted higher albedo). In contrast, in our 
single population hypothesis the lowest mass in colliders is 
« 0.2 M0, and so the Kuiper Belt again appears to be an 
outlier, or at the extreme end of the range. 

4.2 Detectability of Kuiper Belt analogues 

An excess of order one per cent above the stellar photosphere 
will be very difficult to detect around nearby stars, even with 
the newly-launched Herschel mission. The PACS/70 chop- 
ping mode was predicted to detect « 3 mjy sources with 
5(7 confidence in 1 hour, whereas the Sun if seen for exam- 
ple at 6 pc would have a photospheric signal of 135 mJy. 

^ We include only the classical belt and neglect dust that could 
be gene rated over a wider be lt by collisions among scattered disc 
objects. iLevison et al.l ||2008| ) suggest such an event created a re- 
cently discovered family of comets, but models are not well con- 
strained by this one data point. 

^ This value based on the Sun's effective temperature was 
checke d ag ainst a photospheric calculation for the MIPS/70 pass- 
band flVillin g ct al. 2008) for the 'Solar twin' star HD 146233 
dSoubiran fc Triaudll2004i'l . 



Such an excess, at our Solar upper limit of 0.02, could thus 
in principle be robustly detected in an hour. In practice, 
there are only a dozen Sun-like stars within this distance, 
while beyond it a disc of Solar System dimensions would 
span less than three 5-arcsec telescope beams, and so un- 
certainties in subtracting the blended stellar signal would 
be re-introduced. However, Herschel will be able to explore 
somewhat larger excesses very efficiently. For example, if 
5 % uncertainties in stellar models allow the detection of 
i?70 f» 0.15, Herschel is predicted to detect debris in 35 % 
of systems (according to the grey line extrapolation in Fig- 
ure 1), compared to about 20 % of the stars observed with 
Spitzer. 



5 IMPLICATIONS 

The single population hypothesis discussed above fits both 
the debris systems detected by Spitzer and the net upper 
limit for the remaining stars, once the switch of dust removal 
to Foynting-Robertson drag is taken into account. If the Sun 
is a member of this hypothetical population, it must be one 
of the least dusty systems, within the last ten per cent. The 
estimates of mass in colliding comets suggest the Kuiper 
Belt population may actually be a depopulated outlier, with 
only a few-hundredths of an Earth-mass in colliding bodies 
up to tens of kilometre sizes. Roughly a third of similar stars 
are predicted to host an order of magnitude more comets; 
therefore, implications for impacts on any terrestrial planets 
need to be considered. 

Since the Sun is a normal mid-main sequence G-dwarf, 
these results are surprising. Further, the primordial Kuiper 
Belt should have been much dustier and a bright mem- 
ber of the inferred single population. It is estimated that 
at least 10 of rocky mater ial was needed in ord er for 
the Kuiper Belt bodies to form l|Morbidelli et al]|2003l ). and 
this total belt mass would place the young Sun around 
s ~ 0.1 in Figure 1. However, between the formation stage at 
~ 0.1 Gyr and today, the Kuiper Belt s hould only have de - 
clined about four- fold in dust luminosity (|Lohne et aLlbOO"^ ). 
so its present-day rank would be a; ~ 0.3, inconsistent with 
the observed rank of a; > 0.9. 

The key to these results is the dispersal of many Kuiper 
Belt bodi es, so that t he population is now very depleted 
(|Morbidelli et al.ll200l e.g.). The Nice model (|Gomes et all 
I2005D proposes that dispersal occurred when Jupiter and 
Saturn crossed their mutual 2:1 mean motion resonance 
as they migrated outwards; this event would have moved 
Uranus and Neptune onto enlarged eccentric orbits and con- 
sequently disrupted the comet belt. This model can explain 
ejection of many objects, and can be synchronized with the 
Late Heavy Bombardment of the Earth-Moon system at 
about 0.7 Gyr. However, the existence and timing of this 
event d epend on the precise ar chitecture of the system of gas 
giants ilh ommes et al. I [200^ 1. Such as dispersal of around 
99 % of the comet mass could readily shift the Solar Sys- 
tem's excess to a low state and thus a high ran k, as inferred 
here - a process recently modelled in detail bv I Booth et al.l 
(|2009h . However, they find that < 12 % of Sun-like stars can 
have undergone such an event, as there is no global drop in 
dustiness at mature ages (Figure 2), so syste m clearing can 
not have been common. iGaspar et al.l l|2009l ) infer roughly 
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similar proportions, based on the small number of stars that 
could be undergoing clearout events in the 760-Myr old clus- 
ter Praesepe. Thus in the present context, although comet- 
clearing by giant planets could shift the ranks of individual 
stars to much higher values, the data imply that such stars 
are only a small sub-set. This may be neglected in the sin- 
gle population hypothesis, as the net excess at a; > 0.3 is 
presently poorly defined, and so switching a few per cent of 
planet-hosting stars to a low- or zero-dust state would have 
little effect on the net debris estimate. 

5.1 Relation to giant planets 

From long-term tren ds in Doppler wobble surveys, 
ICumming et al.l l|2008h estimate that « 18 % of Sun-hke 
stars host giant planets, of above about Saturn's mass and 
orbiting within 20 AU. Here we assume that any such planet 
orbiting within 3 AU would be disruptive, for example per- 
turbing a terrestrial planet at ~ 1 AU into an unstable or 
eccentric orbit, or in some migr ation scenarios preve nting 
a habitable planet from forming (|Fogg 8z NelsonlboO Q). Ex- 
cluding this 8 % of systems leaves 10 % of stars hosting a 
gas giant at 3-20 AU, among which four out of every ten 
plane ts (generally abo ve Jupiter's mass) are already discov- 
ered (jCumming et al.l[2008 ). A small correction should be 
made for approximately 1 % of star systems found to host 
another lower-mass giant that is closer in, although this is 
less than the uncertainty of about 3 % in the extrapolation 
to orbits as distant as 20 AU. The result is that « 9 ±3 % of 
Sun-like stars should host gas giants no closer in than 3 AU, 
and are therefore reasonably analogous to the Solar Sys- 
tem, with one-third of this population already discovered. 
Hypothesizing optimistically that all these stars could also 
have an unknown ice giant planet like Uranus or Neptune, 
the minimum conditions are met for gas giant interactions 
that could end with disruption of a comet belt outside the 
ice giant (should such a belt exist). 

It is presently uncertain how giant planet systems and 
comet belts are related. Few star s have both gas giants an d 
debris discs ((Greaves et al...2007i : iMoro-Martm et al.ll2007h . 
and systems with both phenomena h ave similar excesses to 
non-planet-hosts (jBrvden et al.|[2009[ ). However, there could 
be a class of systems with many infalling comets and fre- 
quent impacts on any terrestrial planet (s). One archetype 
is HD 10647 which has a very high Ero of 50 and so an 
infer red > 100 Mm o f colliders; this object is > 1 Gyr 
old iLiseau et aLllioOsI ). while the four dusty planet-hosts 
of Figure 3 are about mid-main sequence. Consequently, 
it could be very difficult for complex life to ever evolve 
in these systems with at least twenty times more collid- 
ers than the Sun, if gas giants perturb these bodies into 
the inner regions. (Energetic impacts more frequently than 
ever y ~ 10 Myr would probably hinder recovery of biodiver- 
sity (jKrug fc Patzkowskvll2004l ): however, creation of large 
heated impact ba sins could favour the appea rance of simple 
thermophilic life l|Abramov fc Moizsisll2009l ).) On the other 
hand, whether there is any connection between faint debris 
levels and gas giants is very difficult to determine. There 
are 14 planets hosts in our sample that do not have debris 
detections, and for these the mean E70 is 0.05±0.05, not sig- 
nificantly different from 0.03 ±0.02 estimated for null-debris 
stars in general. 



5.2 Analogue exo-Earths 

An analogue planet to our own is hypothesized to have a 
roughly similar impact rate. Impacts much faster than the 
rate experienced on Earth, of about one every 100 Myr for 
10 km bodies, are likely to cause a serious extinction prob- 
lem (e.g. halving the number of species after each event). 
However, a much slower impact rate could mean that life 
would stay at a very simple level, with no need to adapt to 
changes in environment. Here we consider a cometary popu- 
lation comparable to or lower than that of the Sun as similar, 
but also note that the slope of the inferred E70 is shallow 
at large x, so a considerable population could be only a few 
times dustier (depending on the poorly known Solar -E70). 
With this basic hypothesis, we then need to link impact 
rates to the configuration of the giant planets and comet 
belts. Unfortunately, thi s is computationally e xpensive to 
study for many systems (iHorner fc Jones! [2OO8I . e.g.). Mod- 
elling (Greaves, Jeffers, Horner et al., in prep.) will compare 
the Solar System, with rather few comets but several per- 
turbing giant planets, to representative common extrasolar 
architectures, such as many more comets but fewer perturb- 
ing planets. In the interim, we can consider two cases at 
opposite extremes: one where debris and giant planets are 
unconnected phenomena, and one where gas giants beyond 
3 AU are a hypothetical signpost to systems where comets 
are cleared out. 

The results above then suggest that < 10 % of stars 
have as few comets as the Sun (based on its E70 upper limit) , 
while about 9 % of stars host an innermost gas giant beyond 
3 AU. If these phenomena are uncorrelated, then the product 
of these rates implies that only up to ~ 1 % of Sun-like 
star systems would be analogous to our own, in the sense 
of a small comet population concurrent with outer-system 
giant planets. In the alternative case where giant planets 
can help to enforce a low comet population, this fraction of 
stars rises to an upper limit of 9 %, assuming further-out ice 
giants are also present and so interactions can clear out the 
comet belt. This estimate of 1-9 % for an analogue system 
may be optimistic, if the exo-Earth impact rate is critically 
sensitive to the configuration of gas giants; for example, the 
clearing of the Kuiper Belt required Jupiter a nd Saturn to 
cross a particular resonance, in the model of iGomes et al.l 
|2005l ). Further, there may not be an exo-Earth around all 
Sun- like stars, or not within the Habitable Zone where water 
can be liquid on the surface; see lRavmond et "all (|2007l ) for a 
model prediction that stars of > O.8M0 are the most suitable 
hosts. On the other hand, our estimate of only a few per cent 
of analogue systems is too pessimistic if in fact the single 
population hypothesis does not apply, and there are S> 10 % 
of Sun-like stars with little dust. This is not possible to test 
with the current data, but Herschel observations should soon 
be able to test whether many stars host little dust. 

Adopting up to 9 % as the preliminary estimate of the 
frequency of analogue systems allows us to estimate the dis- 
tance to the nearest counterpart. In conjunction with the 
space density of Sun-like hosts, which for mid-F to late-K 
dwarfs is 0.01 pc~^, the nearest analogue system should lie 
within 13 pc. The nearest presently known analogue candi- 
date (in the sense of planet and comet content) is the G8 star 
HD 154345 at 18 pc distance, which has no obvious Spitzer 
70 [im excess and hosts a Jupiter- like planet at 4 AU. How- 
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ever, if only about one-third of the systems with gas giants 
beyond 3 AU have so far been discovered, a completeness 
correction suggests that a nearer analogue could lie within 
about 12 pc, or slightly closer given that not all nearby Sun- 
analogues as yet have published debris data. 



6 CONCLUSIONS 

Systems roughly analogous to our own, hosting a modest 
comet population and gas giants at a few AU or beyond, 
are inferred to be rather uncommon. If Sun-like stars of- 
ten host terrestrial planets in the habitable zone, only a few 
per cent of these are likely to have a similar debris environ- 
ment. A few systems may have catastrophic bombardment, 
even at mid-main sequence age, if gas giants perturb some of 
their numerous comets into the inner system. At the other 
extreme, the majority population should be stars without 
giant planets but with many comets, for which further mod- 
elling is needed to assess the rate of planetary impacts. The 
closest presently-known system that is roughly like our own 
lies at 18 pc, but another probably exists with around 10 pc 
given the completeness of debris and giant planet surveys 
so far. This result is encouraging for future facilities aiming 
to study habitable exo-Earths, such as DARWIN, TPF and 
ELT. 



ACKNOWLEDGMENTS 

We thank the referee for comments that helped to clarify 
several lines of argument. JSG thanks STFC and SUPA for 
support of this work. 



REFERENCES 

Abramov O., Mojzsis S.J., 2009, Nature 459, 419 
Aumann H.H., Good J.C., 1990, ApJ 350, 408 
Backman D.E., Dasgupta A., Stencel R.E., 1995, ApJ 450, 
L35 

Beichman C.A. et al., 2006a, ApJ 652, 1674 

Bernstein G.M., Trilling D.E., Allen R.L., Brown M.E., 

Holman M., Malhotra R., 2004, AJ 128, 1364 
Booth M., Wyatt M.C., MorbideUi A., Moro-Martm A., 

Levison H.F., 2009, arXiv:090£3755] 
Bryden G. et al., 2009, ApJ 705, 1226 
Bryden G. et al., 2006, ApJ 636, 1098 
Carpenter J.M. et al., 2009, ApJS 181, 197 
Gumming A., Butler R.R, Marcy G.W., Vogt S.S., Wright 

J.T., Fischer D.A., 2008, RASP 120, 531 
Decin G., Dominik C, Malfait K., Mayor M., Waelkens C, 

2000, A&A 357, 533 
Dominik C, Decin G., 2003, ApJ 598, 626 
Fogg M.J., Nelson R.R, 2009, A&A 498, 575 
Gaspar A. et al., 2009, ApJ 697, 1578 

Gomes R., Levison H.F., Tsiganis K., MorbideUi A., Nature 
435, 466 

Greaves J.S., Wyatt M.C., Bryden G., 2009, MNRAS 397, 
757 

Greaves J.S., Fischer D.A., Wyatt M.C., Beichman C.A., 
Bryden G., 2007, MNRAS 378, LI 



Greaves J. S., Wyatt M. C, Holland W. S., Dent W. R. 

F., 2004, MNRAS 351, L54 
Greaves J.S., Wyatt M.C., 2003, MNRAS 345, 1212 
Habing H.J. et al., 2001, A&A 365, 545 
Hahn J.M., Malhotra R., 2005, AJ 130, 2392 
Hillenbrand L.A. et al., 2008, ApJ 677, 630 
Horner J., Jones B.W., 2008, A&G 49, 1.22 
Jeffers S.V., Mauley S.P., Bailey M.E., Asher D.J., 2001, 

MNRAS 327, 126 
Kim S. et al., 2007, BAAS 39, 814 

Krug A.Z., Patzkowsky M.E., 2007, Paleobiology 33, 435 
Landgraf M., Liou J-C, Zook H.A., Griin E., 2002, AJ, 
123, 2857 

Levison H.F. MorbideUi A., Vokrouhlicky D.D., Bottke 

W.F. 2008, AJ 136, 1079 
Liseau R. et al., 2008, A&A 480, L47 

Lohne T., Krivov A.V., Rodmann J., 2008, ApJ 673, 1123 

Meyer M.R. et al., 2008, ApJ 674, L181 

Meyer M.R., Backman D.E., Weinberger A. J., Wyatt M.C., 
2007, in Protostars and Planets V, B. Reipurth, D. Jewitt, 
and K. Keil (eds.). University of Arizona Press, Tucson, 
p.573 

MorbideUi A., Brown M.E., Levison H.F., 2003, Earth, 

Moon & Planets, 92, 1 
Moro-Martm A. et al. 2007, ApJ 658, 1312 
Raymond S.N., Scalo J., Meadows V.S., 2007, ApJ 669, 606 
Soubiran C, Triaud A., 2004, A&A 418, 1089 
Stapelfeldt K.R. et al., 2004, ApJS 154, 458 
Su K.Y.L. et al, 2006, ApJ 653, 675 

Thommes E.W., Bryden G, Wu Y., Rasio F.A., 2008, ApJ 
675, 1538 

Trilling D.E., Beichman C.A., Bryden G., Rieke G.H., Su 

K.Y.L. , 2008, ApJ 674, 1086 
Trilling D.E. et al., 2007, ApJ 658, 1289 
Wyatt M.C., 2008, ARAA 46, 339 
Wyatt M.C., 2005, A& 433, 1007 

Wyatt M.C., Smith R., Su K.Y.L., Rieke, G.H., Greaves 

J.S., Beichman C.A., Bryden G., 2007, ApJ 663, 365 
Wyatt M.C., Dent W.R.F., 2002, MNRAS 334, 589 
Zahnle K., Arndt N., CockeU C, Halliday A., Nisbet E., 
Selsis F, Sleep N.H., 2007, Space Science Reviews 129, 35 

This paper has been typeset from a L^jX/ M^^jK file prepared 
by the author. 



